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Abstract 

In this paper we review the theory of the "falhng chimney", which deals with the breaking in 
mid-air of tall structures, when they fall to the ground. We show that these ruptures can be caused 
by either shear forces, typically developing near the base, or by the bending of the structure, which 
is caused primarily by the internal bending moment. 

In the latter case the breaking is more likely to occur between one third and one half of the height 
of the chimney. Small scale toy models are used to reproduce the dynamics of the falling chimney. 
By examining photos taken during the fall of these models we test the adequacy of the outlined 
theory. This type of experiment, easy to perform and conceptually challenging, can become part 
of a rotational mechanics lab for undergraduate students. 
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I. INTRODUCTION 



One of the most interesting demonstrations for an introductory mechanics course is the 
"Falhng Chimney - Free Fall Paradox," as it was named by Sutton in his classical book 
Demonstration Experiments in Physics} In the original version of this demonstration a ball 
is placed at one end of a uniform stick, which is pivoted at the other end and makes initially 
an angle of about 30° with the horizontal. The elevated end of the stick is suddenly dropped, 
together with the ball, thus showing a very counter-intuitive behavior. The falling end of 
the stick accelerates at a greater rate than the free-falling ball, proving that its acceleration 
is greater than g, the acceleration of gravity. 

A simplified version of the experiment can be performed just with a meter stick and a 
coin. The stick is supported in the horizontal position by two fingers, placed near the two 
ends. A coin is set on the stick near one end, which is suddenly released. The effect is similar 
to the previous demonstration: the falling end of the rotating stick eventually acquires an 
acceleration greater than that of the freely falling coin, which loses contact with the stick 
surface and lags behind the falling stick. 

A description of the first version of the experiment can be found in almost every book of 
physics demonstrations, ^'^'^ sometimes with different names ("Falling Stick", "Hinged Stick 
and Falling Ball", and others). Photographic descriptions or even video-clips of this demo 
can be found on-line in several web-pages (see our web-page^S for a collection of related 
links) . 

In addition, countless papers exist in the literature; we have traced several of these, from 
the 1930's to the present. Some of the earliest discussions can be found in Constantinides- 
and Ludekrf (as well in the book by Sutton^), followed by many others . ^i^i^^i^-^i^^ These 
concentrate mostly on the simple explanation of the effect, which relies on the concept 
of "center of percussion" of the rotating stick (a simple introduction to this concept can 
be found in Bloomfielcii^). This particular point of the stick (located at a distance from 
the hinged end equal to two thirds of the length, for a uniform stick) is moving with the 
same acceleration as a particle under gravity, constrained to move along the same circular 
path. Points on the stick beyond the center of percussion descend with accelerations greater 
than that of particles freely moving under gravity, on their respective circular paths. As a 
consequence of this, if the initial angle formed by the stick with the horizontal is less than 
about 35°, the end point will possess at all times a vertical component of the acceleration 
greater than g, producing the effect described above. 

Several variations of the basic demonstration also exist , ^^i^^i'^^i'^'^'i'^^i'^^ the majority of which 
suggest attaching an additional mass to the rotating stick at different positions. The effect 
for the student or the viewer is even less intuitive than the original version: an additional 
mass placed near the end of the stick actually reduces the acceleration of the end point, 
affecting substantially the outcome of the experiment. In general the addition of a mass at 
any point on the stick will increase both the total torque on the system (thus increasing the 
rotational acceleration) and the moment of inertia of the system around the axis of rotation 
(resulting in a decreased rotational acceleration). The center of percussion of the stick still 
plays a key role: if the additional mass is placed beyond it, the effect of the increased moment 
of inertia dominates and the acceleration of the rotational motion will be reduced. If the 
mass is placed before the center of percussion, the increase in the torque will dominate and 
the rotational motion will be enhanced. The effect is null if the mass is placed exactly at 
the center of percussion (a complete discussion of this effect can be found in Bartletfe^" and 
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Haber-Schaimi^) . 

The next logical step is to analyze the behavior of a real falling chimney. Almost in- 
variably a tall chimney, falling to the ground like the stick in the previous discussion, will 
break in mid- air at some characteristic height. This is well documented in several photos 
reproduced in the literature, such as the one which appeared on the cover of the September 
1976 issue of The Physics Teacher (other photos can be found in Bundy^'' and Bartlett,— 
or also on our web-page^). The causes of such breaking, the height of the rupture point 
and the angle at which the breaking is most likely to occur, are the most natural questions 
which arise. 

The first analysis^ of this problem compared the fall of the real chimney to the fall of 
the hinged stick, but wrongly identified the center of percussion (at about two thirds of the 
height) as the probable point of rupture. Reynolda^^ first identified the possible causes of 
the breaking with the shear forces and the bending moment originating within the structure 
of the toppling chimney. More detailed analyses were given by Bundy^S and Madsen24 
(the most complete papers we found on the subject) while simplified explanations are also 
reported . ^^i^^i^^ It even appears in graduate student study guides ,-^^ although the chimney 
is shown bending the wrong way in one of these books. 

In this paper we review the theory of the real falling chimney, outlined by Madsen,— 
aiming for a complete and clear explanation of this phenomenon in Sects. ITTland inTl Then, 
in Sect. IIVI we propose simple ways of using small scale models (literally toy models - made 
with toy blocks and bricks) to test effectively the outlined theory. More information on these 
toy models can also be found on our web-site,— together with photos and movie clips of the 
experiments we have performed. 

II. ROTATIONAL MOTION OF THE FALLING CHIMNEY 

The rotational motion of a falling chimney under gravity is equivalent to that of the falling 
hinged stick of Sect. HJ We can describe it as in Fig. Q where we use polar coordinates r and 
6 (with and as unit vectors) for the position of an arbitrary point A on the longitudinal 
axis of the chimney, measuring the angle 6 from the vertical direction. We treat the chimney 
as a uniform rigid body of mass m and height H, under the action of its weight W = mg, 
applied to the center of gravity (basically the center of mass -CM- of the body), and a force 
F exerted by the ground on the base of the chimney, assumed to act at a single point (we 
neglect air resistance, or any other applied force). In plane polar coordinates: 



The moment of inertia of the chimney can be approximated with the one for a uniform 
thin rod, with rotation axis perpendicular to the length and passing through one end:— 



Applying the torque equation 19 = r^, for a rotation around the origin, with an external 
torque given by = mg^ sin 6*, we find the angular acceleration 



W = WrCr + Weee = —mg cos 9er + mg sin 9ee 
F = FrCr + Feee- 



(la) 
(lb) 



/ = -mH^. 
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FIG. 1: The falling chimney described as a rotating uniform stick. The external forces are the 
weight of the body applied to the center of gravity, and the constraint force at the base. 

A simple integration, using = ^ = %0 = f^^^^^' S^^^^ ^^e angular velocity 

^' = 3|:(l-cos0), (4) 

assuming that the chimney starts moving from rest and is initially in the vertical direction. 
A further integration of Eq. 01 can lead to 9{t) in terms of elliptic integrals. 

■ 2 ^ 

We recall that the acceleration in polar coordinates can be written as a=r=(r — r9 )er + 
{r9 + 2r9)ee, so that, for a point yl at a fixed distance r from the origin, it becomes 

a = arCr + agee = —r9 + r9eg. (5) 

For a point at two thirds of the height, r = combining Eqs. El and El we get ag{r = 

= ^H9 = gsm9, proving that this particular point is the center of percussion of the 
body, as already mentioned in Sect. HI 

The torque equation allowed us to determine the angular acceleration of the motion in 
Eq. El We can use this result and Newton's second law for the motion of the center of mass 
(CM) of the whole chimney to determine the unknown force F at the base. The vectorial 
equation is 

mrcM = W + F (6) 
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FIG. 2: The forces acting on the lower portion of the chimney, due to the upper part and the 
action of the constraint at the base, are shown here. The two insets explain the definition of the 
bending moment in terms of a couple of forces. The resulting deformation of the structure is also 
shown for the two possible cases. 



which, for r = y, sphts into radial and angular equations, 

H -2 

—m—9 =Fr — mgcos9 (7a) 
H 

m—9 = Fg + mg sin 9. (7b) 
Using Eqs. El and IH the two components of the force F are easily determined: 

Ft = ^mg (cos 9 — | ) (8a) 



2 V 5^ 
Fq = —-mg sin 9. (8b) 



III. INTERNAL FORCES AND BENDING MOMENT 



We now move to the analysis of the internal forces which develop inside the structure 
of the falling chimney. The resulting stresses and bending moment are the causes of the 
rupture of the toppling chimney. Consider, as in Fig. |21 an arbitrary lower portion of the 
chimney of height r (as opposed to the total height H) and the forces acting on this part of 
the structure due to the upper portion and the base. The weight of the lower portion is now 
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W{r) = mgjj (assuming again a uniform structure, so that the weight is proportional to 
the height of the considered portion) and it is applied to the center of gravity of this lower 
portion (at a distance | from the origin). In polar coordinates: 

W(r) = Wr{r)er + Wg{r)e0 = —mg— cosOcr + mg— sin^eg. (9) 

H H 

The force F, applied at the base, is still the same as in Eqs. IHil EEl but we have to add 
the action of the upper part on the lower portion. We follow here the general analysis of 
the internal forces and moments which can be found in every textbook on Statics (see for 
exampl o ^ '^ ) and which can be easily adapted to our case. 

The distribution of all the internal forces, at the cross section being considered, can be 
equivalently described by a resultant force and a resultant moment acting at a specific point 
of the cross section (typically the "centroid" of the sectioned area, in our case simply the 
central point of the section, on the longitudinal axis). In particular, the resultant force can 
be decomposed into a transverse shearing force S = Seee, and a longitudinal stress force 
(tension or compression) P = PrCr, applied as in Fig. 12 at the cross section between the 
upper and lower portions, and assumed positive in the eg, direction respectively. 

In addition, we have to consider the resultant moment of the forces at the cross section, 
which is usually called the "bending moment" N^, because its effect will ultimately result in 
bending the structure. It is represented in the picture by the curved arrow. Since we treat 
this as a plane problem, the bending moment can only have a component perpendicular to 
the plane of the figure, i.e., in the z direction. No other components are considered here, in 
particular we assume that no torsional moment exist in the structure, which would tend to 
twist the chimney around its longitudinal axis. 

The bending moment N;, = A^^e^ can be thought as originating from a couple of forces, 
f and — f , that can be regarded as applied to the leading and trailing edge of the structure, 
at the cross section considered. This couple of forces is shown explicitly in the papers 
by Bundy2S and Madsen,^^ but we prefer to use directly N;, in our treatment, because 
the bending moment is the result of the whole distribution of forces at the cross section 
considered. The two small insets inside Fig. |21 explain the definition of the bending moment 
in terms of the couple of forces f and — f. We also show the resulting deformation of 
the structure due to a "clockwise" (diagram a), or a "counter-clockwise" bending moment 
(diagram b). The latter case will be the actual deformation of the falling chimney. Nh will 
be assumed to be positive if it acts as in the figure, i.e., a positive component of the torque 
in the z direction (we assume here the use of a right-handed system of coordinate axis). In 
the following we will refer to A^^^ as the bending moment, acting on the lower portion of the 
chimney. '^^ 

Again, we will consider the torque equation and the second law for the motion of the 
center of mass (located at |) for just the lower portion of the chimney (of mass jjItl). It is 
better to analyze the CM motion first. The vector equation 

m^icM = W(r) + F + P + S, (10) 
si 

will split into the radial and angular directions, 

7^2. 2 3 ^ 5 / 3\ 

— m 6 = — mgfl — cos^) — - = —niq — cos^H — mg ( cos^ I + -Pr (Ha) 

2H 2 ' 2 Sy ^ ^ 

m—r:0 = -mqiAnO-— = mq—smO — -mqsmO + Se, (lib) 
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having used Eqs. El IH Eil EEl andlHl We can solve for the longitudinal and transverse forces 



Pr = '1^9 {l-^) [(5 + 3-^)cos^-3(l + -^)] (12a) 

3 / 4 r l\ 

S-, = -mf^sin^ —--- + - , (12b) 



4 ^ 3H ' 3^ 

which depend on the fraction of height jj, the angle of rotation 9, and also the total weight 
mg. Following the analysis by Bundy^^ we plot these two forces in Figs. 3 and 4 respectively, 
normalized to the total weight mg, as a function of the height fraction, for several angles. 

From Fig. El we see that Pr is negative (a compression) for smaller angles, but eventually 
becomes positive (a tension) for angles greater than about 45°. Pr also depends critically 
on jj (for ^ = 0°, Pr represents simply the compression due to the weight of the upper part 
acting on the lower part). This longitudinal force will be combined later with the bending 
moment to determine the total stress at the leading and trailing edges, which is the most 
typical cause of the rupture. 

In Fig. IHwe plot the (transverse) shear force Sg, which can be the other leading cause 
of rupture. It is easily seen that, for any considered angle, the magnitude of the shear force 
1 5*61 1 has an absolute maximum at = (and a positive value), meaning that large shear 
forces, in the eg direction, usually originate near the base. The shear force is always zero 
at one third of the height, and \Sg\ also has a (relative) maximum at |if (with a negative 
value, therefore Sg is in the —'eg direction), but this value is smaller than the one near the 
base. 

From this analysis, it is typically concluded that if the structure breaks from shear stress 
alone, this is usually more likely to happen near the base. This can be seen for example in 
the already mentioned cover photo of the September 1976 issue of The Physics Teacher,^^'^ 
showing the fall of a chimney in Detroit. The two ruptures at the bottom are likely due 
to shear forces, while a third rupture can be seen at about r = OATH, and this is due to 
the combination of bending moment and longitudinal force Pr, as we will explain in the 
following. More photos and detailed pictorial descriptions of chimney ruptures can be found 
in the paper by Bundy.^^^ 

The "bending moment" A^^ can be calculated from the torque equation I{r)9 = r^, where 

now /(r) = ^^j^r'^ is the moment of inertia of just the lower part. 9 will come from Eq. El 
and the total external torque is now = ^Wg{r) + rSg + N^. Using Eqs. lUland lT^Bl we can 
solve the torque equation for the bending moment, obtaining 

Nh = —-mg sm9 r [1 — — ) , (13) 
4 \ H J 

or, in a non-dimensional form 

* 4-^^^l-4)^ (14) 



mgH 4 H \ H 

which is plotted in Fig. 5, as a function of the height fraction and for various angles. 

Nb is always negative, showing that it is actually directed in the opposite way of Fig. El 
(or as in diagram b of Fig. |21). This particular direction of the bending moment will induce 
a tension in the leading edge of the chimney and a compression in the trailing edge. The 
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FIG. 3: The longitudinal stress foree per unit weight of the chimney is shown as a function of the 
height fraction for several angles. Positive values indicate tensions, while negative values represent 
compressions. 
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FIG. 4: The transverse shear force per unit weight of the chimney is shown as a function of the 
height fraction for several angles. Positive values are for forces in the eg direction. 
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FIG. 5: The bending moment, divided by the weight and the height of the chimney, is shown as a 
function of the height fraction for several angles. 
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structure will bend accordingly, with the concavity on the side of the trailing edge, and will 
eventually break in the way shown by the many existing photos. We can see that, for any 
angle, the bending moment is obviously zero at the bottom and at the top of the chimney, 
while it assumes its maximum absolute value at exactly one third of the height H. The 
bending moment alone would therefore induce a rupture at one third of the structure, but 
the total longitudinal stress at the leading edge is also due to the force Pr, as we will show 
next. 

Another interesting relation between the bending moment A^;, and the shear force Sg is 
that they are in general related by a simple derivative, i.e., Se = as it is easy to check 

from Eq. I12bl and Eq. ^1^^ This is a well known relationship of the Statics of beams and 
other structural members (for a complete proof see for example Hibbeler^°). It is a direct 
consequence of the equilibrium equations applied to an infinitesimal longitudinal portion of 
the beam: the change in bending moment along the beam is always equal to the shear force 
applied to that portion of the beam. 

Finally, we can combine A*";, and Pr to compute the total longitudinal stress on the cross 
sectional area between the lower and upper parts. We follow the theory of elasticity and 
deformations in beams, which can be found in classical treatises such as Sommerfelci^ and 
Landau-Lifshitz,— or again in Bundy's paper.— The longitudinal stress is maximum at 
the leading and trailing edges, located at the maximum distance from the longitudinal 
(centroidal) axis of the chimney which lies within the "neutral surface" of the structure, the 
surface which is neither stretched nor compressed. 

For simplicity, we will only consider from here on, structures with uniform square cross 
section of side a, as this is the case of the toy models described in Sect. IIVI In this case the 
stresses at the leading and trailing edge, ai and a-r respectively, can be evaluated from 

Pr aNb , . 

^L/T = J^T^ (15) 

(the upper sign is for ctl, the lower for ax) where A = a? is the area of the square cross 
section of side a, with the factor | representing the distance between the longitudinal axis, 

4 

considered as the "neutral axis," and the two edges. J = ^'is the moment of inertia of the 
cross sectional area computed about the neutral axis (see Sommerfelctfi^^ for details). Using 
also the expressions for Pr and A^f,, we obtain 



(7L/T _ _\(^_ ^ 



mg 2 V H 



5 + 3— j cos^-3(^l + — 



i^^sin^^fl-^)', (16) 



where we normahzed aL/r, dividing by in order to obtain a dimensionless quantity which 
is plotted in Fig. El as a function of the height fraction and for several angles. 

This quantity depends also on the ratio which for a real chimney is of the order 
^ > 10. For the toy models described in Sect. IIVI the value of this ratio is even bigger: 

~ 24 — 61, enhancing the contribution of the second term of Eq. which comes from 
the bending moment Nb- In Fig. IHlwe show the plot for ^ = 24, but similar figures can be 
obtained for different values of the ratio. 

In Eqs. IT31ITB1 and in Fig. 10 the total stresses are considered positive if they represent 
tensions, negative if they are compressions. It is easily seen from the figure that the stress at 
the leading edge 0"^, is initially a compression, but eventually becomes a tension, constantly 
increasing for larger angles; ax on the contrary is usually a compression. It is therefore the 
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FIG. 6: The normalized longitudinal stress at the leading edge (solid) and at the trailing edge 
(dotted), are shown as a function of the height fraction and for several angles {H/a = 24). The 
maxima of the leading edge stress curves are marked by solid points. 
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combination of these intense tensile stresses in the leading part of the chimney (and also 
compressions on the trailing side) that causes the rupture of the chimney. 

This type of breaking is more likely to occur at the positive maximum value of ai- This 
maximum value depends critically on the jj ratio, for a certain angle of rupture 6*, so it 
is possible, just by looking at the maxima of the solid curves in Fig. IHl (marked by solid 
points), to roughly match the height of the rupture point to the angle at which the breaking 
started to occur. As far as the actual prediction of the point of rupture, this would obviously 
depend on the building materials and the construction of the chimney or tower, an analysis 
of which goes beyond the scope of this work. 

It is interesting to note, from Fig. IHl again, that the stress is not always maximum at 
one third of the height (as for the bending moment of Fig. EJ. For small angles of about 
^ ~ 5° — 20° it reaches a maximum for ^ 0.4 — 0.5, while for larger angles it approaches 
the typical ratio ~ |. This means that if the chimney breaks early in its fall, for small 
angles, it is more likely to break near the center; on the contrary if the rupture occurs at 
larger angles, the breaking point is usually shifted toward one third of the height. This is 
the effect of the Pr term in Eq. E[ which modifies the position of the maximum (or 0"^). 

As already mentioned, another factor to be considered is the ratio ^, which is usually in 
the range ^ ~ 5 — 20 for real chimneys, but can be increased up to about ~ 100 with 
our toy models. In Fig. [7| we plot the maxima of the leading edge stress curves, for several 
values of the ratio - = 5, 10, 20, 30, 100 (for the group of curves between the values of 30 
and 100, the parameter is increased by 10 units for each curve). 

These maxima are plotted as continuous curves, showing the corresponding values of the 
ratio jj and the angle 6*, at which the structure is more likely to break, for a given value 
of ^. In other words, connecting for example the solid points of Fig. we would obtain a 
corresponding continuous curve in Fig. [7| with the angle of probable rupture on the vertical 
axis, instead of the normalized stress. 

The dependence of the rupture point on these quantities can be noted in several of the 
existing photos of falling chimneys, when the breaking occurs due to the bending of the 
structure and not for the transverse shear stress near the base mentioned at the beginning 
of this section. The angle of rupture can be roughly estimated by measuring the angle the 
upper part of the chimney forms with the vertical in the photos. This angle tends in fact 
not to change much after the rupture, since the upper part falls without much additional 
rotation. For example the photos in the paper by Bartletlj^'^i refer to chimneys with ^ ^ 10, 
breaking at about 9 ~ 20° — 25° and for ~ 0.47, consistent with the values of Fig. [71 for 
the ^ = 10 curve. Similar behavior can be seen in other photos, since real chimneys tend 
to break very early in their fall, due to the intense stresses originating within their structure. 

IV. TOY MODELS 

In this section we discuss our efforts to reproduce the effects described above, with the 
help of toy models of the falling chimney. These models were constructed with simple toy 
blocks of different type and size, and their fall was filmed with a digital camera, so that we 
could analyze the events frame by frame, to test the theory. Complete details on the type of 
blocks used, experimental settings and video-capture techniques, as well as the complete set 
of our video-recordings and still photos can be found on our web-site,^ and in an upcoming 
publication.— 
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FIG. 7: The maxima of the leading edge stress curves are shown as continuous functions, for several 
values of the ratio H/a. They represent the points, in terms of the height ratio r/H and angle 0, 
at which the structure is more likely to break, due to bending. 
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Wa = 24 
rH =0.50 
theta<10 deg. 



FIG. 8: The first toy model made with wooden blocks. The structure appears to break at = 0.5, 
and at a small angle 9 < 10°, which can be estimated by measuring the angle formed by the upper 
part of the tower with the vertical. 

Bundy noted in his work^ that the use of a model to test a real chimney would be 
useless due to a "scale effect." The stresses inside the chimney depend roughly on the scale 
of the object, so that real chimneys would develop bigger stresses than equivalent small- 
scale models, therefore breaking earlier in their fall. Nevertheless we found it interesting to 
reproduce these effects in small scale models to test especially the discussion based on Figs. 
iniand[7| It was also supposed to be difficullj2& to show these effects with toy models.— 

Fig. IHlis the first example of one of our toy models. We made a tower by simply stacking 
24 wooden toy blocks of cubic shape, for a total height H = 0.76 m, mass m = 0.32 kg 
and a ratio ^ = 24, the value used in Fig. [Ul The tower was set into the falling motion 
by removing a support at the bottom, inducing a rotation without slipping at the bottom 
point. The picture clearly shows the "rupture" due to bending of the structure at exactly 
half the height, jj = 0.50, and for a small angle 6 < 10°, which again can be estimated by 
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FIG. 9: The second toy model made with plastic blocks. The structure appears to break at 
■jj = 0.40, and at a larger angle 6 ~ 30° — 35°, which can be estimated by measuring the angle 
formed by the upper part of the tower with the vertical. 

measuring the angle the upper part of the chimney forms with the vertical direction. This 
is in good agreement with the position of the maximum for the solid-10° curve in Fig. Inland 
also with the data of Fig. [7J using the ^ = 20 curve, showing that the theory is applicable 
also to these small-scale models. 

Our second example is a taller tower {H = 1.9 m, m = 0.65 kg, and ^ = 61) made with 
100 plastic blocks of a very popular brand of toy bricks. The blocks are inserted on top of 
each other so that bending of the structure is allowed, but shear stress cannot possibly break 
the tower. The 100 toy blocks are arranged by color to subdivide the structure into three 
equal parts, and also the position of the center is marked. This time the rupture occurs for 
an angle around 9 ~ 30° — 35°, and at the height ratio jj = 0.40. This is consistent with the 
30° solid curve in Fig. IHl while the data from Fig. [7| (for ^ = 60) would suggest a smaller 
angle of rupture. It is actually difficult, with this type of toy bricks, to estimate the angle 
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at which the structure begins to bend. 

We performed several other experiments, varying the dimensions of the towers, the type 
of blocks, always obtaining results consistent with the theory. We can conclude that it is 
actually easy to reproduce the bending and breaking of chimneys with small scale towers, 
and this type of experiment could be made part of an undergraduate laboratory class for 
rotational mechanics, with some minor adaptations and changes.— 

V. CONCLUSION 

In this paper we reviewed the theory of the falling chimney, showing that the rupture can 
be caused by either shear forces typically near the base, or by the bending of the structure 
which is caused primarily by the bending moment, but is also affected by the longitudinal 
stress force. In the latter case the breaking is more likely to occur between one third and 
one half of the height of the chimney. 

This point of rupture is also related to the angle of rupture and this relationship can be 
verified in the many existing photographic reports of falling chimneys. We also constructed 
several small scale toy models, showing that it is possible with them to reproduce the 
dynamics of the fall. By examining photos taken during the fall of these models we were 
able to confirm the theoretical model outlined in this paper. 
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